Multi-band and nonlinear hopping corrections to the 3D Bose-Fermi-Hubbard model 
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Recent experiments revealed the importance of higher-band effects for the Mott insulator (MI) - superfluid 
transition (SF) of ultracold bosonic atoms or mixtures of bosons and fermions in deep optical lattices [Best et 
al, PRL 102, 030408 (2009); Will et al, Nature 465, 197 (2010)]. In the present work, we derive an effective 
lowest-band Hamiltonian in 3D that generalizes the standard Bose-Fermi Hubbard model taking these effects as 
well as nonlinear corrections of the tunneling amplitudes mediated by interspecies interactions into account. It 
is shown that a correct description of the lattice states in terms of the bare-lattice Wannier functions rather than 
approximations such as harmonic oscillator states is essential. In contrast to self-consistent approaches based on 
effective Wannier functions our approach captures the observed reduction of the superfluid phase for repulsive 
interspecies interactions. 

PACS numbers: 03.75.Lm, 03.75.Mn, 37.10.Jk, 67.85.Pq 



I. INTRODUCTION 

Ultracold atoms in optical lattices provide unique and 
highly controlable realizations of various many-body Hamil- 
tonians [1-6]. Theoretical descriptions of these systems in 
the case of deep lattice potentials usually employ lowest-band 
models only [1,7]. However, it was found recently that for lat- 
tice bosons with strong interaction contributions to the Hamil- 
tonian beyond the single-band approximation with nearest- 
neighbor hopping and local two-particle interactions need to 
be taken into account [8]. E.g., using the method of quan- 
tum phase diffusion, the value of the two-body interaction 
U for bosons in a deep optical lattice was measured directly 
and found to deviate from the prediction of the tight-binding 
model derived in [1]. These experiments also revealed the 
presence of additional local three- and four-body interactions 
not accounted for in the single-band Bose-Hubbard Hamil- 
tonian. A perturbative derivation of these terms based on 
harmonic-oscillator approximations was given by Johnson et 
al. [9]. 

In the case of boson-fermion mixtures, the situation is more 
involved. The first experiments on mixtures with attractive 
interspecies interaction [10, 11] displayed a decrease of the 
bosonic superfluidity in the presence of fermions. This initi- 
ated a controversial discussion about the nature of the effect. 
Explanations ranged from localization effects of bosons in- 
duced by fermions [11, 12] to heating due to the admixture 
[10, 13]. Numerical results also predicted the opposite be- 
havior, i.e., the enhancement of bosonic superfluidity due to 
fermions [14] with a more detailed discussion in [15]. The sit- 
uation remained unclear until a systematic experimental study 
of the dependence of the shift in the bosonic SF - MI transi- 
tion on the boson-fermion interaction [16] and the subsequent 
observation of higher-order interactions in the mixture. This 
shows, that again higher-order band effects need to be taken 
into account. 

The influence of higher Bloch bands in the Bose-Fermi 
mixture can be described by two different approaches: In 
the first approach one assumes that the single-particle Wan- 
nier functions are altered due to the modification of the lattice 
potential for one species by the interspecies interaction with 



the other [17], which is then calculated in a self-consistent 
manner. The agreement of these results to experimentally 
observed shifts of the SF-MI transition is very good for the 
case of attractive boson-fermion interaction (see [16]). The 
method fails however for repulsive interactions where exper- 
iments showed contrary to intuition again a reduction of su- 
perfluidity [16]. Besides this shortcoming, the self-consistent 
potential approach has a conceptual weakness as it can only 
be applied close to the Mott-insulating phase. The second 
approach to include higher bands is an elimination scheme 
leading to an effective single-band Hamiltonian similar to the 
pure bosonic case [9, 18, 19]. This approach, although techni- 
cally more involved, is more satisfactory from a fundamental 
point of view. It did not result in quantitatively satisfactory 
predictions so far, however. We will show here that this is 
because (i) an important non-linear correction to the hopping 
mediated by the inter-species interaction and present already 
in absence of higher-band corrections has been missed out and 
(ii) harmonic oscillator approximations to the Wannier func- 
tions which have been used before, lead to gross errors when 
considering higher band effects. 

We here present an adiabatic elimination scheme for Bose- 
Fermi mixtures obtained independently from [9, 18, 19], re- 
sulting in an effective first-band BFH-Hamiltonian [20]. In 
contrast to [9] and [18, 19] we use correct Wannier functions, 
which will be shown to be essential. Furthermore we find that 
already within the lowest Bloch band the inter-species interac- 
tion leads to important nonlinear corrections to the tunneling 
matrix elements of bosons and fermions. For a fixed number 
of fermions per site, the effective Hamiltonian is equivalent 
to the Bose-Hubbard model with renormalized parameters U 
and J for which expressions are given in a closed form. This 
allows for a direct study of the influence of the boson-fermion 
interactions on the bosonic superfluid to Mott-insulator tran- 
sition within this level of approximation. It is shown that non- 
linear hopping together with higher-band corrections lead to a 
reduction of the bosonic superfluidity when adding fermions 
for both, attractive and repulsive inter-species interactions. 

The outline of the present work is as follows. After deriv- 
ing the general multi-band Hamiltonian of interacting spin- 
polarized fermions and bosons in a deep lattice in the follow- 
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ing section, we introduce the first important addition to the 
standard BFHM in section III, the nonlinear hopping correc- 
tion. Restricting to leading contributions, we derive an effec- 
tive single-band Hamiltonian by adiabatic elimination of the 
higher bands in section IV. Finally, using the resulting gener- 
alized BFHM, the effect of a varying boson-fermion interac- 
tion is studied in detail in section V. 



H. MODEL 

In 3D, ultracold Bose-Fermi mixtures in an external poten- 
tial are described by the continuous Hamiltonian [7] 
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where the index b (f) at the field operators i& refers to bosonic 
(fermionic) quantities and V b (r) [V^(r)] is the external po- 
tential consisting of possible trapping potentials as well as 
the optical lattice V£ t (r) = ?] b X>in 2 (fcar Q ) [ViL( r ) = 
rjf ^2 sin 2 (k a r a )]. The intra- and interspecies interaction 
constants are defined as 
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with ma = being the reduced mass and at>b/bf the 

intra- and interspecies s-wave scattering length, respectively. 

Whereas in the standard approach the field operators in (1) 
are expanded in terms of Wannier functions for the first band 
only, we here use an expansion to all Bloch bands: 



*/W=EE^*-j)- 



(3) 



The operator b v - } denotes the annihilation of a boson 

(fermion) in the i/-th band at site j and Wu ? (r — j) is the 
corresponding Wannier function of the v-th band located at 
site j. The vector v = {VxiVy, v z} denotes the band index. 
The Wannier functions factorize as 



w b Jf(v)=w b J/(x) w b J/(y)w b J/(z) 
whith the one-dimensional Wannier function w b J* (x) 



(4) 



Using the expansion of the field operator, the full multi- 
band Bose-Fermi-Hubbard Hamiltonian can be expressed as: 

jl-j4 
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The generalized hopping amplitudes (still containing local en- 
ergy contributions) 
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and the generalized interaction amplitudes 
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(9) 



are defined as usual. In the following we restrict our model 
in such a way, that only the most relevant terms are kept. 
Note that many of the matrix elements vanish because of 
the symmetry of the Wannier functions [21]. Unless stated 
otherwise we restrict ourselves to local contributions in 
interaction terms, i.e. ji = ■ • ■ = j 4 in U^J 4 and V^-J 1 
and in this case we drop the site indices. 

With these restrictions, the general multi-band Hamiltonian 
can be cast in the following form 



Ujtl 



where the first term 



Hi — i^BFHM + -ffnlii 



(10) 



(11) 



describes the (pure) first-band (1 = {1,1,1}) dynamics 
consisting of the standard Bose-Fermi-Hubbard part ffeFHM 
[7] and nonlinear hopping corrections i? n im which will be 
discussed in the next section. The second term H® incorpo- 
rates the (free) dynamics within the v-th band and H u ^ grT 
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describes the coupling between arbitrary bands v>, fi, g, a. 
The prime in the sum indicates that at least one multi-index 
has to be different from the others. This general form of the 
full Hamiltonian serves as the starting point of our study. 



III. NONLINEAR HOPPING CORRECTION 

Even when virtual transitions to higher bands are disre- 
garded there are important corrections to the standard BFHM 
if the boson-fermion interaction V becomes large. The in- 
terspecies interaction term in (1) gives rize to a correction to 
the bosonic (and fermionic) tunneling amplitude proportional 
to the occupation number of the corresponding complemen- 
tary species. These contributions, in the following termed as 
nonlinear hopping contributions, have been considered before 
[22, 23], but have been missed out in earlier discussions of 
corrections to the BFHM [18, 19]. 

To establish notation let us recall first the usual single-band 
BFHM 

&=-j£&JSj + |5>j(A|-l) 

<ij> v (12) 

<ij> " J 

The amplitudes are determined by 

TT — TTlM') T/ — T/jjjJ t = ji+^-i 7= _ 7J + S 'j 
u — U 111H v — 'UlU J — J l,l j J — J l,l 

with e being an unit vector in one of the three lattice direc- 
tions. Due to the isotropic setup, the choice of the direction is 
irelevant. From eq. (8) and (9) two types of nonlinear hopping 
corrections arise: From the boson-boson interaction we obtain 

<ij> 

whereas the boson-fermion interaction leads to both, bosonic 
and fermionic hopping corrections: 

<ij> „ (14) 
+ fi(ni + nj)fr 

<ij> 

The corresponding nonlinear hopping amplitudes read 

jb — rrj+ej.jj jf _ T/j+e.jJ j J T/jJ,j+e,j 
J nl — '-'llll ' J nl — r llll ' J nl — * 1111 

Since we are interested in the influence of the fermions to 
the bosons we assume in the following the fermions to be ho- 
mogenously distributed. This assumption also used in [16, 17] 
proved to be valid in the trap center and gives a consider- 
able simplification. This amounts to replacing the fermionic 
number-operators by the fermionic filling: mj — > m. Further- 
more, the bosonic density-operators in eqns. (13) and (14) are 



replaced by the filling of the Mott-lobe under consideration, 
fij — > n, for simplicity. 

Alltogether, this allows us to write a Hamiltonian including 
corrections from the nonlinear hopping contributions. Defin- 
ing the effective bosonic hopping amplitude as 

J[n, m] = J — 2n — fnJja, (15) 

the system is recast in the form of a pure BHM with density 
dependend hopping: 

H ctf = -J[n, m] E «ta + (16) 

<ij> J 

Analyzing the resulting predictions for the MI-SF transition 
as a function of the filling and the interspecies interaction (see 
figure 3) one recognizes a substantial reduction of bosonic su- 
perfluidity for increasing interaction on the attractive side and 
a corresponding enhancement on the repulsive side, showing 
the importance of nonlinear hopping terms for the precise de- 
termination of the MI-SF transition. Compared to the exper- 
imental results [16], two main points arise. First, although 
pointing into the right direction for attractive interactions, the 
overall shift is too small compared to the experimental ob- 
servation. Second, for repulsive interactions, the transition is 
shifted to larger lattice depths, in contrast to the experimental 
findings. 

IV. EFFECTIVE SINGLE-BAND HAMILTONIAN 

In the following we derive an effective single-band Hamil- 
tonian that takes into account the coupling to higher bands. 
The derivation is structured in the following way: We use an 
adiabatic elimination scheme presented in appendix A which 
reduces the main task to the calculation of the the second 
order cumulant ((T Hj(t + T)Hi(t))) in the interaction 
picture, where the average is taken over the higher bands. 
The full interaction Hamiltonian Hi = Yl' H ullQcr is 

then reduced according to the relevant contributions of the 
cumulant. Finally, a reduction of the effective bosonic 
scattering matrix (A5) gives the full effective single-band 
Bose-Fermi-Hubbard model. 

When calculating the cumulant ((T #/(r + T)Hi(t))) in 
(A5), the interaction Hamiltonian of the full multi-band Bose- 
Fermi-Hubbard model can be reduced considerable. Keeping 
only terms that lead to non-zero contributions in lowest or- 
der, it is easy to see that only those terms in Hj matter, where 
particles are transfered to higher bands by -ffr(r) and down 
again by Hj(t + T). In the following we restrict ourselves to 
precisely those contributions and furthermore treat only local 
contributions since these are dominant. Three relevant pro- 
cesses are found: 

1 . Single particle transitions to a certain band v 
{1,1,1,1} o {v, 1,1,1} 
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FIG. 1 . (Color online) Matrix elements for the coupling of the higher 
Bloch bands to the first band via the generalized interaction (9). 
Bosonic contributions from (8) are equivalent. Bosons are shown 
as orange circles and fermions in black. V u tui describes the tran- 
sition of a boson and a fermion from the first (higher) to the higher 
(first) band; V^xvX gives two particles (boson and fermion) which 
perform a transition to bands v and fi. V^m derscribes a fermion- 
mediated single particle transition of a boson, where Vn^,i is a 
boson-mediated transition of a fermion. 

These contributions can be understood as density- 
mediated band transitions, where the matrix elements 
U vlll , V^ X11 , V lwll are only non-zero for odd bands 
v. [24] Note that from now on, the upper site-indices 
are omitted if they are all the same. 

2. Double-transition to the same band v 

{1,1,1,1} O {!/,!/, 1,1} 

In this situations, two particles undergo a transition to 
the same band and all bands are incorporated. The ma- 
trix elements are U vv \\ and V vv \\. 

3. Double-transition to different bands v, /j, 



{1,1,1,1} o K/M,!} 

In this combined process, the two different bands have 
to be both either even or odd with matrix elements 
Uvpn and V vl ixi 

The remaining important contributions to the full multi- 
band BFHM result from the kinetic energy of the particles. 
Restricting to the usual nearest neigbour hoppings within a 
given Bloch band (v = fi and |ji — J2I = 1) and the energy 
of the particles within a band (v = fi and ji = j2), these are 

4. the band energies A^, and A/, 

5. the intraband nearest-neighbor hopping for bosons J„ 
and correspondingly for the fermions J v . 

Hopping between sites with |ji — ja| 7^ 1 is omitted since it 
is unimportant. In appendix B, the different contributions to 
the Hamiltonian as well as the hoppings and band energies 
are defined in detail. Figure 1 gives a sketch of the different 
contributions taken into account. Shown are only processes 
involving fermions. 

From the effective bosonic scattering matrix in (A5), the 
effective single-band BFHM is derived by applying a Markov 
approximation [25]. This amounts to replacing first-band op- 
erators at time r + T by the corresponding operators at time r 
which is valid since the timescale of the higher-band dynam- 
ics is much shorter than in the first band because of the larger 
hopping amplitude [26]. The resulting Hamiltonian is lengthy 
and shows the full form is given in appendix C. 

The effective Hamiltonian (CI) contains non-local interac- 
tion and long-range tunneling terms. These result from virtual 
transitions into higher bands and subsequent tunneling pro- 
cesses in these bands. As these terms rapidly decrease with 
increasing distance |d| between the involved lattice sites, it is 
sufficient to take into account only the leading order contribu- 
tions, i.e. only local interaction terms (|d| =0) and only near- 
est neigbour hopping (|d| = ±1). This leads to the following 
extensions compared to the standard single-band BFHM: 



- ff f ^3 ^3 ^2 V2 \ 



(17) 



]T Afrj + Aim, J i< hj,m, mj] h £ ft J[n u ftj] f] + £ { J™ (btfb* + j( 2 ) SJfl/jS,} 

j J (ij> <ij> <ij> L J 



Here some new terms arize, for instance correlated two- 
particle tunneling J^> and J^ 2 \ Most prominent is the 
appearance of the three-body interactions Us and Vij.The 
bosonic has recently been measured by means of quantum 
phase diffusion [8]. It should be noted that in the experiments 
in [8] also higher order nonlinear interactions were detected. 



Since our approach is only second order in the interaction- 
induced intra-band coupling, these terms cannot be repro- 
duced however. Beside the new terms, the higher bands lead to 
a renormalization of the usual single-band BFHM parameters. 
Whereas the local two-body interaction amplitudes U2 and V2 
only depend on the band structure, the hopping amplitudes are 
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altered, leading to density mediated hopping processes. For 
the bosonic ones, the hopping now is of the form 

J[hi, hj,rhi, toj] = J - J nl (nj + Hi) 

2^ (™j + ^0 + a 

+ /3 rhihj + 7 hirhj + <5 rhimj (18) 

and the density dependence is directly seen. For all parameters 
occuring in (17), full expressions can be found in appendix D. 



V. INFLUENCE OF FERMIONS ON THE BOSONIC MI-SF 
TRANSITION 

In order to discuss the phase transition of the bosonic sub- 
system, we make further approximations. Coming from the 
Mott insulator side of the phase transition, the local num- 
ber of bosons is approximately given by the integer average 
filling, i.e., (nj) sw n. For the fermionic species, we also 
replace the number operator by the average fermion number 
raj — >• m = 1, assuming a homogeneous filling of fermions 
in the lattice. Having an experimental realization with cold 
atoms in mind, this is a valid assumption in the center of the 
harmonic trap at least for attractive inter-species interactions. 
It should be valid however also for slight inter-species repul- 
sion. This assumption is also supported by the results of [16], 
where the actual fermionic density did not influence the tran- 
sition from a Mott-insulator to a superfluid (for medium and 
large filling). It also agrees with the result in [17] which is 
based on this assumption, and which shows a good agree- 
ment to the experimental results. All further contributions 
in the Hamiltonian such as the bosonic three-particle interac- 
tion and two-particle hoppings are neglected in the following. 
With these approximations, the renormalized Bose-Hubbard 
Hamiltonian for the n-th Mott lobe with mean fermionic fill- 
ing to can be written as 

= -J[n, m] £ b\k + ^ E "J (*1 - !) < 19 > 

(ij> " J 

with 

J[n, m] = J- 2n J h nl - m J f nl (20) 

x — ^ * / Tfl \ ^ 

- 2^ T b.» { U vlll H + V ulll If J 

U[m] = U 2 +mV 3 (21) 

The final form of the bosonic Hamiltonian will now be used 
to discuss the influence of the boson-fermion interaction on 
the Mott-insulator to superfluid transition. Following the ex- 
perimental procedure presented in [16], we consider the shift 
of the bosonic transition as a function of the boson-fermion 
interaction determined by the scattering length cibf, with a 
special emphasis on repulsive interaction where no theoreti- 
cal prediction exists so far. 



nonlinear bosonic 

nl bosonic + bands 

Ml to SF transition 

: / / 


✓^harmonic approximation/ 
/ pure bosonic / i^'' 


nl bosons +^ands + fermions-^/' 


termionic shift 



10^ — 1 1 ^ ^ 1 ' — 1 " 1 1 1 

5 6 7 8 9 10 11 12 13 14 15 



potential depth ti b [E R( J 

FIG. 2. (Color online) Ratio of effective interaction U to effective 
tunneling rate J for unity fermion filling m — 1 and Mott lobe with 
n — 1 as function of normalized lattice depth rft for the bosons, 
and for attractive boson-fermion interaction with a scattering length 
dbf = —400 ao- The horizontal dotted line gives the critical value 
for the MF - SF transition point in the Bose-Hubbard model. Shown 
are the harmonic oscillator approximation together with different lev- 
els of corrections as described in the main text based on exact Wan- 
nier functions. Assumed is a perfect match between fermionic and 
bosonic Wannier functions, r\F = t]b- 



To determine the transition point, we calculate the bosonic 
hopping (20) and interaction amplitude (21) using numeri- 
cally determined Wannier functions. The knowlegde of the 
critical ratio U / J of the MI to SF transition from analytic 
or numerical results [27-29] allows for the precise localiza- 
tion of the transition point [30]. This method is displayed in 
figure 2 where the ratio of the effective interaction strength 
U[\] and the effective tunneling rate J[3, 1] as per eq. (20) to 
(21) are plotted as a function of the normalized lattice depth 
■qb, which describes the amplitude of the periodic lattice po- 
tential of the bosons V^ t in units of the recoil energy of the 
bosons £^ ec = ft 2 k 2 / (2m;,). As indicated, unity fermion fill- 
ing to = 1 is assumed and the bosonic Mott lobe with n = 3 is 
considered. The horizontal dotted line gives the critical value 
for the MI - SF transition [27] and the crossing of this line 
with the different curves, which illustrate the relative contri- 
bution of the various correction terms, determines the poten- 
tial depth at which the phase transition occurs. The different 
levels of approximation shown in figure 2 are 

1. harmonic oscillator: 

plain BHM, harmonic oscillator approximation 

2. pure bosonic: 

plain BHM, proper Wannier functions 

3. nonlinear bosonic: 

BHM extended by nonlinear (bosonic) hopping correc- 
tion 

4. nonlinear bosonic with higher bands: 

inclusion of all bands with v a < 25; this gives the ref- 
erence point for the shift of the transition 
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5. nonlinear bosonic and fermionic with higher bands: 
inclusion of fermions; nonlinear hopping correction and 
higher bands (y a < 25) 

One clearly recognizes a substantial shift of the transition 
point to lower potential depth in qualitative agreement with 
the experiment. It is also apparent that using harmonic os- 
cillator approximations leads to a large error of the predicted 
transition point. This shows that the use of the correct Wan- 
nier functions is crucial for obtaining reliable predictions. 

Figure 3 shows the shift of the MI - SF transition point for 
the first four lobes as a function of the boson-fermion scat- 
tering length abf. The solid lines include all corrections de- 
scribed earlier, where the amount of the shift is measured rela- 
tive to the nonlinear bosonic case including higher bands, i.e., 
relative to the real bosonic transition point. Thus the figure 
corresponds to the shift of the transition point when fermions 
with unity filling are added to the system. For each Mott lobe 
three curves are shown corresponding to different ratios of 
Vf/ 7 lb which illustrates the effect of different masses and/or 
different polarizabilities of the bosonic and fermionic species 
as discussed in Appendix E. The dashed-dotted curves give 
the contributions of the (first band) nonlinear hopping correc- 
tions only (bosons and fermions). One recognizes that for in- 
creasingly attractive interactions between the species there is 
an increasing shift of the transition point towards smaller po- 
tential depth, corresponding to a reduction of bosonic super- 
fluidity in the presence of fermions. Interestingly one recog- 
nizes that for repulsive interspecies interactions, virtual tran- 
sitions to higher Bloch bands tend to counteract the effect of 
the fermion induced nonlinear tunneling. For larger values 
of n there is again a shift of the MI - SF transition point to- 
wards smaller lattice depth, i.e. again a reduction of bosonic 
superfluidity! The latter effect has both been observed in the 
experiments [16], but has not been fully understood so far. 
In the calculations, the bands are summed up to a maximal 
multi-index !/ max = {25, 25, 25}, including altogether 15625 
bands. For this number of bands, a satisfactory convergence 
of the effective amplitudes U and J is found. Overall, our 
second order approach inlcuding the nonlinear corrections al- 
ready provides an intuitive explanation for the behaviour of 
the system in the experiment. This especially holds for the re- 
pulsive case, where the agreement to the experimental results 
is on a quantitative level. 



VI. SUMMARY AND OUTLOOK 

In the present paper we studied the influence of nonlinear 
tunneling processes and higher Bloch bands on the dynam- 
ics of a mixture of bosons and fermions in a deep optical lat- 
tice in a full 3D setup. Taking into account virtual inter-band 
transitions in lowest non-vanishing order and contributions of 
the originally continuous interaction to tunneling processes 
we derived an an effective lowest-band Hamiltonian extend- 
ing the standard Bose-Fermi Hubbard model. This Hamil- 
tonian contains interaction-mediated nonlinear corrections to 
the tunneling rates, renomalized two-body interactions, and 
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FIG. 3. (Color online) Shift of the bosonic Mott-insulator to super- 
fluid transition as a function of the boson-fermion scattering length 
at, f for different Mott lobes (solid lines, n = 1 . . . 4, from bottom to 
top) in one dimension. The gray-shaded region depicts the influence 
of a mismatch of the bosonic and fermionic lattice depth. The dot- 
dashed lines give the shifts of the transition solely from the nonlin- 
ear tunneling corrections. Dashed horizontal lines give the transition 
points for the pure bosonic system. 



effective three-body interaction terms. We showed that an ac- 
curate determination of the effective model parameter requires 
the use of the correct Wannier functions of the correspond- 
ing single-particle model. As differences in the tails of the 
wavefunctions are essential, the use of approximate harmonic 
oscillator wavefunctions can lead to large errors. The effec- 
tive model allows for a study of the effect of admixing spin- 
polarized fermionic atoms to the bosonic superfluid to Mott- 
insulator transition when changing the boson-fermion inter- 
action strength. Our model recovers qualitatively all features 
observed in the experiment. In particular we found that bo- 
son superfluidity is reduced both for attractive and repulsive 
inter-species interactions. The latter has not been reproduced 
so far with other methods such as the self-consistent potential 
approach. 

It should be noted that our model does not take into account 
heating effects and effects such as phase separation due to the 
presence of an inhomogeneous trapping potential, which have 
recently been shown to significantly affect the MI-SF transi- 
tion point already in the lowest band [31, 32]. We thus expect 
that a complete picture of the experimental observations will 
require a proper inclusion of higher-band effects and nonlin- 
ear tunneling as derived in the present paper, as well as effects 
from heating and a trapping potential. Finally it should be 
mentioned that our approach is limited to the second order in 
intra-band processes. In higher-order perturbation theory ef- 
fective four-body, five-body, etc. interactions will arise, which 
play however a less and less important role. Nevertheless, we 
expect that the higher orders should substantially improve the 
results, especially for repulsive interactions. 
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APPENDIX 

Appendix A: Adiabatic elimination scheme 

As long as the interaction energies U and V as well as 
the temperature are small compared to the band gap between 
lowest and first excited Bloch band, the population of higher 
bands can be neglected. However, as noted before, there are 
virtual transitions to higher bands which need to be taken into 
account. In the following we employ an adiabatic elimina- 
tion scheme of higher Bloch bands starting from the general 
multiband Hamiltonian (10). This scheme, which is also used 
in [33] for the Bose-Fermi-Hubbard model in the ultrafast- 
fermion limit, is equivalent to degenerate perturbation theory 
[27, 34] and allows for a proper description of the reduced 
system. For this, the Hamiltonian (10) is split up into a free 
and an interaction part H = Hf ree + Hi with 



Hf ree = H i + 



ujtl 



H®, 



(Al) 



Hi = h ^,q,t- (A2) 

Transforming to the interaction picture, the dynamics of the 
free part is incorporated by the time dependent interaction 
Hamiltonian Hi(t) = e~i H i r " T Hi f rBeT . Adiabatic 
elimination is carried out for the time evolution operator (scat- 
tering matrix) of the full system given by 



S = T cxp \ - - 



&tHi{t) 



(A3) 



We now trace out the higher-band degrees of freedom, assum- 
ing empty higher bands. Using Kubo's cumulant expansion 
[35] 



S 

<exp{,sX}) x = exp<[ ((X m )) 



k m—1 



(A4) 



up to second order in the interband coupling, the effective 
scattering matrix for the lowest band reads 



S oS = Tcxpi 



(A5) 



As discussed in section IV, the different terms to the Hamil- 
tonian are given by 

1 . Single particle transitions to a certain band v 



+U llul b\b v b{bi + U XXXv b{b v b{bi- 

+v ulll bikfih + v lvll b\bjih- 

+V llvl blkftfi + V lllM blkflfl 

2. Double-transition to the same band v 

2 J2 [ U wi.i fytyMh + U lluv b{b{b u b v + 
j 

+V vXvl bthftfi + Vx v xv b[bJlU 

3. Double-transition to different bands is, /x 

jEKmii t>lblb 1 b 1 + U 1 lvn b[b\b v bfj,+ 
j 



(Bl) 



(B2) 



(B3) 



Only local contributions are taken into account and thus the 
spatial index j is ommited for the moment. The further intra- 
band contributions are defined as 



4. the band energy 



K= j d 3 r<(r) 



2m, 



-A + V x (r) 



<(r), (B4) 



5. the intraband nearest-neighbor hopping for bosons 

h 2 



J v = J d 3 rwl(r-e) 



2mi, 



A + V\v) 



<(r) (B5) 



and correspondingly for the fermions J v . 



AT((THi(t + T)Hi{t))) 



The first order does not lead to any contributions because of 
the vacuum in the higher bands and due to the nature of the 
interband couplings. Obviously the effective bosonic Hamil- 
tonian is connected to the second order cumulants of operators 
in higher Bloch bands, ((A B)) = (A B) — (A)(B) [35]. 



Appendix C: Full effective flrst-band BFHM 

Under the assumptions made in the main text (i. e., only 
local contributions, nearest-neighbour hopping, etc.), the fi- 
nal form of the effective Hamiltonian is found from equa- 
tions (A5) together with the interband couplings from (Bl) 
to (B3)in Markov approximation. This yields 
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Hf = H x 



^1 jd I 



^i/Hl^i/111^6,1/ 
2 



«j+d /jVd^j/j f 



E eI%^<4U\ 2 



jd 



I (jWO T d ft ft Lf.l 

(CI) 

In the Hamiltonian, the time integrals over the bosonic and 
fermionic correlators are defined as 



T d 



h Jo 



dT(b^ +d (T + T)bl.(r)), (C2) 



x (/w+dCr + T)/^^)), (C3) 

and correspondingly and lf h . The two-point correla- 
tion functions of bosons and fermions in the iv-th band read 



(S I /, j+d (r + r)S+ J (r)>, = -l^ 



e - 2 «Ii e i r 4'" (C4) 



e- 2 ^ e ^ T <\ (C5) 



Carrying out the time integration gives in the thermody- 
namic limit, which is obtained for L —> oo by setting £ = ¥ 
and changing jij^ k to JJJ d 3 £ yields: 



I d = -L_ 

x '" (2tt) 3 



1 + S VU S, 



Here 



(2tt 



|0 



d 3 £ d 3 £' 



-i^-dp— tf'-d 



(C6) 



6^(0 + 6^^')' 

(C7) 



e t,I '(f)= E 2J„ a cos(&) + A> a (C8) 

cn—x,y,z 

e f '"(0= E 2X Q cos(C a )+A{ Q (C9) 



is the energy of a boson respectively fermion in the higher 
band and x distinguishes between bosons (x = b) and 
fermions (x = /). 



Appendix D: Definition of constants in Hamiftonian (17) 



As used in Hamiltonian (17), the full expressions of the 
different parameters are: 



Density-mediated fermionic or bosonic hopping: 



J[h h n i+& ,m h m i+& ] = 



J. 



J - Jul (Aj+8 + Aj) - ~f (™j+e + m- 3 ) 



E x t{owr 

f#l I 



^J/lll^J/111 - - . (^I/lll) - ~ 1 /r-.is. 

o n j+z m i + - A m i+& m- 2 ). (Dl) 



J, 



E 

^i 



nj+gnj, (D2) 



pair tunneling amplitude: 



U- 



j+e j+eJJ 
1111 



E 



will 



(D3) 



+ E 



vfill; bb.ufji 



J(2) 



T/j+e,j,j+e,j (T/ ^2 J' 

•'llll , l V i/li/l7 x ( 



E 

\p (V^i M i) 2 T § 



(D4) 
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renormalized two-particle interactions: 

(D5) 

v 2 = v + J2 ^ f ' vv 




(D6) 



E 



^ -w 



bf,ufii 



three-body interactions 

(7 3 = 6E(^iii) 2 ^. 



ujtl 



(D7) 
(D8) 



Appendix E: Lattice effects 

The lattice potentials for bosons and fermions are both cre- 
ated by the same laser field and the only externally control- 
lable parameter is the intensity of this lattice laser. In or- 
der to see how the parameters of the effective lattice model, 
such as tunneling rates and interaction constants depend on 
this laser intensity one needs to take into account that there 
is always a fixed ratio / between the bosonic and fermonic 
potential depths for given atomic species and transitions. To 
determine / we note that the optical lattice is generated by an 
off-resonant standing laser field. The potential itself results 
from the ac-Stark shift. As shown in [36], it is given by 



Vpot(r) 



3^c 2 



2T, 



J(r) (El) 



J 0,D 1 



J 0.Z5 2 



in rotating wave approximation for a typical alkali D-line dou- 
blet, where each line contributes independently if the laser is 
sufficiently far detuned from the atomic transitions. The im- 
portant parameters are the decay rates Td 1 2 of the excited 
states, Ad 1 2 = o;i ascr — ujq.Di 2 me detunings of the laser fre- 
quency wiascr from the atomic transition frequencies ujq,Di 2 
and I(r) = Iq sin 2 (kr) the laser intensity. 

Conveniently, all energies in the system are normalized to 
the recoil energy of the bosonic species given by E% ec = y-^-. 
The wavenumber k depends on the chosen optical lattice. 
The (normalized) lattice potential for the bosons thus reads 



V[at( r ) = t]b sin 2 (kr). It is useful to rewrite the optical lattice 
potential for the fermionic atoms with respect to the bosonic 
optical lattice as V^Jr) = r]f sin 2 (kr), where r]f = f rjb- 
From eq. (El) we find 



21' 



(E2) 



At this point, we specify the experimental system. In the 
previous discussions, we analyzed the experiment reported in 
[16] and use the parameters given there. A mixture of bosonic 
87 Rb and fermionic 40 K is cooled and put into an optical lat- 
tice with ojj = 755 nm. For Rubidium and Potassium, the 
transition wavelengths and decay rates are given by 



766.5 nm 



= 38.7 x 10 6 Hz 



Rb 
r Rb 



'D 2 



769.9 nm 



Rb 
D 2 



795.0 nm 
36.1 x 10 6 Hz 
780.2 nm 



(E3) 



r£ 2 = 38.2 x 10 6 Hz 



a 



T%\ = 38.1 x 10 6 Hz. 



Using these values, / in equation (E2) evaluates to / = 
2.04043, which means, that the fermionic lattice potential, 
in terms of the bosonic recoil energy is twice as deep as the 
bosonic one. For the calculation of the Wannier functions 
of bosons and fermions one has to take into account how- 
ever also the different masses of the particles. Expressing the 
Schrodinger equation for the single-particle fermionic wave- 
function $/(r) in terms of the bosonic quantities rjb and m,b, 
one finds 



2m,h 



m f 7 2/i \ 

— fVb sin (kr) 



$,(r) = -±E<f> f (r). 
m b 

(E4) 

One recognizes that the difference between the fermionic 
Wannier functions and the bosonic ones is determined only 
by the factor ^-/. Since for the experiment in [16] 



E b , 



= 2.175 



(E5) 



the factor / is almost compensated, —f-f = 0.93. Thus 
the bosonic and fermionic Wannier functions are to a good 
approximation identical with a maximal overlap. Neverthe- 
less, figure 3 also display results including a mismatch of 
the bosonic and fermionic Wannier functions, depicted by 
the gray shaded regions. The upper (lower) boundary on 
the attractive side and the lower (upper) boundary on the 
repulsive site corresponds to the results for a mismatch of 
^-f = 0.7 C^-f = 1.3), indicating the importance of a 
good control of the mismatch in the precise determination of 
the transition shift. 
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